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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


MEMORANDUM 4-19-59A 


SUPERSONIC AND MOMENT -OF -AREA RULES COMBINED FOR RAPID 
ZERO-LIFT WAVE-DRAG CALCULATIONS 
By Lionel L. Levy, Jr. 


SUMMARY 


The concepts of the supersonic area rule and the moment-of-area rule 
are combined to develop a new method for calculating zero-lift wave drag 
which is amenable to the use of ordinary desk calculators. The total 
zero -lift wave drag of a configuration is calculated by the new method 
as the Siam of the wave drag of each component alone plus the interference 
between components . In calculating the separate contributions each 
component or pair of components is analyzed over the smallest allowable 
length in order to improve the convergence of the series expression for 
the wave drag. The accuracy of the present method is evaluated by com- 
paring the total zero-lift wave-drag solutions for several simplified 
configurations obtained by the present method with solutions given by 
slender-body and linearized theory. The accuracy and computational time 
required by the present method are also evaluated relative to the super- 
sonic area rule and the moment-of-area rule. 

The results of the evaluation indicate that total zero-lift wave- 
drag solutions for simplified configurations can be obtained by the 
present method which differ from solutions given by slender-body and 
linearized theory by less than 6 percent. This accuracy for simplified 
configurations was obtained from only nine terms of the series expression 
for the wave drag as a result of calculating the total zero -lift wave 
drag by parts. For the same number of terms these results represent an 
accuracy greater than that for solutions obtained by either of the two 
methods upon which the present method is based, except in a few isolated 
cases. For the excepted cases, solutions by the present method and the 
supersonic area rule are identical. Solutions by the present method are 
obtained in one fifth the computing time required by the supersonic area 
rule. This difference in computing time of course would be substantially 
reduced if the complete procedures for both methods were programed on 
electronic computing machines. 
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INTRODUCTION 


The supersonic area rule presented in reference 1 provides a useful 
tool for calculating the zero-lift wave drag at supersonic Mach numbers 
of configurations consisting of slender bodies , thin wings, and thin 
tails. The zero-lift wave drag of a given configuration was shown to 
be the average of the wave drag of a series of equivalent bodies of revo- 
lution. Since the development of the supersonic area rule, various 
numerical methods, easily adaptable to punch-card computing machines, 
have been developed for calculating the zero-lift wave drag of these 
equivalent bodies of revolution. To mention a few, references 2, and k, 
respectively, present methods based on a knowledge of the area distribution, 
the first derivative of the area distribution, and the second derivative 
of the area distribution of the equivalent bodies of revolution. 

In reference 5 a method was developed for calculating zero-lift wave 
drag at low supersonic Mach numbers which does not require a knowledge 
of the equivalent bodies of revolution, but uses, instead, a convenient 
set of geometric parameters which are solely a function of the area and 
moment distributions of the given configuration. The geometric parameters 
consist of double moments of the area distribution of the given config- 
uration calculated about both the streamwise and spanwise axes of the 
configuration. For this reason the method of reference 5 is referred to 
as the moment -of -area rule. With this method, it is practical to make 
zero-lift wave-drag calculations with ordinary desk calculators. 

Similar to the supersonic area rule, the moment-of-area rule also 
evaluates a series expression for the zerc-lift wave drag. As a result 
of a simplifying assumption the moment-of- area-rule series expression 
for the wave drag does not converge as raj icily as does that for the super- 
sonic area rule. The accuracy and the number of double moments (terms of 
the series) required by the moment-of-arec rule to obtain adequate con- 
vergence of the series at supersonic Mach numbers are too great for the 
method to be practical with desk calculators. At a Mach number of 1, 
however, the number of moments (calculated about only the spanwise axis 
of the given configuration) required are greatly reduced to the point 
where desk calculators are practical. It was reasoned, therefore, that 
if a convenient and systematic technique could be developed for calculating 
the moments of equivalent bodies of revolt tion, the concepts of the 
moment-of-area rule at Mach number 1 and the supersonic area rule could 
be combined to develop a method for calculating zero-lift wave drag for 
a wide range of Mach numbers by means of crdinary desk calculators. 

It is the purpose of this report to presert the development and evaluation 
of such a method. The method will be evaluated by comparing results 
computed by the method developed herein with those obtained by other 
existing methods of calculation. 



LIST OF IMPORTANT SYMBOLS 


An ( 0 ) 9 ) 

b 

B 

c 

d 


D(P) 

D(M) 


K 


l 

HM) 

M 

Mpk 


M ok<M) 


coefficients of a Fourier sine series expansion of S T (x,[3,0) 
wing span 
1 - A 

wing chord at the vertical plane of symmetry 

lateral distance between the longitudinal axes of a pair of 
bodies of revolution 

zero -lift wave drag of a configuration 

zero-lift wave drag of an equivalent body of revolution of a 
confi guration 


tangent of the sweep angle of the 50 -percent chord line of any 
sheared panel of a wing or tail surface in dimensionless 
coordinates 

length of a body of revolution 

length of an equi valent body of revolution 

free -stream Mach number 

moments of the area distribution of a given configuration 
(see eq. (9)) 

moments of the area distribution of an equivalent body of 
revolution 


N 


number of terms or harmonics used in the calculation of wave 
drag 


q free -stream dynamic pressure 

S(xj(3,0) frontal projection of the area distribution intercepted on a 
given configuration by a set of parallel oblique planes 
tangent to the Mach cones 

t Q wing maximum thickness at the vertical plane of symmetry 

t ( x , y ) thickness distribution of a configuration 
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x,y,s Cartesian coordinates in the free -stream , spanwise, and 

thickness directions , respect j. vely 

a local dimensionless airfoil-section chord station measured 

from the local 50 -percent chord 


P 

8 

G 

A 


T 


n/m 2 - 1 

longitudinal distance "between the lateral axes of a pair of 
sheared bodies of revolution 

angle defining the orientation of the parallel oblique planes 
tangent to the Mach cones (see sketch (a)) 

plan -form taper ratio 

dimensionless Cartesian coordinates in the free-stream, span- 
wise, and thickness directions, respectively 

maximum thickness ratio 


t(cx) dimensionless thickness distribution of the local airfoil 

section 


t ( | , rj,K) dimensionless thickness distribution of one panel of a sheared 

plan form 




spanwise variation of the dimensionless airfoil-section thick- 
ness along lines of constant percent chord 


Subscripts 


k,m,n,\ 
P,r,s J 

I 

L 

T 

U 

W 


indices of summation 

in t e rf e r e n c e c on di t i on s 

lower panel of a plan form, r\ < 0 

total 

upper panel of a plan form, r\ > 0 
wing panel alone conditions 
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properties of an exposed wing measured at or relative to the 
average chord at the wing-body juncture 

E,H component (body) E or component (body) H 

EH,UL interference , specifically, between bodies E and H or upper 

and lower wing panels, respectively 


Superscripts 

differentiation with respect to a coordinate in the free- 
stream direction, except as noted in equation (C2) 

dimensionless symbol 


METHOD 

Basic Methods 


Before the present method for calculating zero -lift wave drag is 
developed, it is well to review briefly the basic concepts of the two 
methods upon which the present analysis is based. 

Supersonic area rule .- In reference 1 the zero-lift wave drag of 
configurations consisting of slender bodies, thin wings, and thin tail 
surfaces was shown to be the average of the wave drag of a series of 
equivalent bodies of revolution. This fact is given analytically by 


D(P) = | 



n=i 


n[A n (p,0)] 2 d0 


which, for convenience, is written 


(i) 


D(0) = i / D(p,0)d0 

2* J Q 


( 2 ) 


where 


d(m) = J q 


CO 

^ n[A n (p,0)] 2 


n=i 


( 3 ) 
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As defined in reference 1 the term D(p,0) is the zero-lift wave drag, 
as the Mach number approaches 1, of each e iui valent body of revolution. 
For a given free-stream Mach number, or P, each value of 0 specifies 
one member of the series of equivalent bodies of revolution. The 
functions A n (3,0) are coefficients of a Fourier sine series expansion 
of the first derivative of the area distribution of the equivalent bodies 
of revolution and are defined by 


A n (P,e) = 


-TC 


dS(x,p ,a : 

5x 


sin(ncp)dcp 




where 


x 


Up, g) 

2 


COS Cp 


( 5 ) 


The normal cross-sectional area distribution of the equivalent bodies of 
revolution, S(x,p,0), is obtained as the frontal projection of the area 
distribution intercepted on the given configuration by a set of parallel 
oblique planes tangent to the Mach cones. The term 2(3,0) is the length 
of the equivalent bodies of revolution. Tie coordinate system, angles, 
and Mach planes are defined in sketch (a) . 



Sketch (a) 

Application of the basic methods is restricted to a particular group 
of configurations. This group satisfies tie following conditions: 
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i. Body components and wing and tail components of a configuration 
are, respectively, sufficiently slender and sufficiently thin that only 
negligible errors are introduced into the calculation of S(x,p,0) by 
assuming the oblique planes to be normal to the horizontal plane. 

ii. The area and first derivative of the area distribution of each 
equivalent body of revolution must not have discontinuities. 

iii. The slope of the area distribution at the ends of each equiva- 
lent body of revolution must be zero. 

The accuracy of the results obtained by the supersonic area rule 
depends upon the accuracy with which the Fourier coefficients are evalu- 
ated and upon the convergence of the series expression for the wave drag 
of each equivalent body given by equation (3). This latter dependency 
results from the fact that in practice the infinite series must be 
terminated at some finite number of terms. In references 2 and 3 the 
Fourier coefficients are evaluated from a knowledge of the area distri- 
bution and the first derivative of the area distribution of the equivalent 
bodies of revolution, respectively. In both references N = 25 provides 
adequate convergence of the series for configurations which satisfy the 
restrictive conditions listed above. In reference 4 it was pointed out 
that for configurations which satisfied the above restrictive conditions 
but also had singularities at the ends of the second derivative of the 
area distribution of at least one component of the equivalent bodies of 
revolution, a larger number of terms was required to provide adequate 
convergence of the series. In this case all of the 49 available terms 
were required. Finally, equation (3) must be evaluated for enough values 
of 0 to define D(P,0) for integration. 

Moment -of -area rule . - Zero-lift wave-drag calculations by the moment- 
of-area rule differ basically from those by the supersonic area rule in 
the manner in which the Fourier coefficients are evaluated. Unless all 
parts of a configuration lie within the nose Mach cone and the forward 
Mach cone from the tail, the equivalent body length, z(3,0), will be 
greater than the actual length of the given configuration for some 
values of 0. As shown in reference 5> however, by considering stream- 
wise body extensions of vanishingly small cross-sectional area, one can 
asstime the length of each equivalent body of revolution to be constant 
and equal to or greater than the length of the longest equivalent body, 
say l . In this manner equation ( 5 ) becomes 

x = z/2 cos cp (6) 

As a result of this assumption the moment -of -are a series expression for 
the zero-lift wave drag does not converge as rapidly as does that for the 
supersonic area rule. Upon substitution of equation (6) in equation (4) the 
Fourier coefficients can be expanded in a finite series expressible in 
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terms of a convenient set of geometric parameters, Mp k , and. powers of 
(3 cos 9 such that 

2 n " 2 

A n (p,e) = 2 (j) Y, cos 0)P U) 

p=o 


where 


n “f 2 -p-* 

^np = ^ Snpk (jJ Mpk 

k=o 


( 8 ) 


These geometric parameters are double moments of the thickness distri- 
bution of the given configuration calculated about both the streamwise 
(x) and spanwise (y) axes of the configuration. The origin of the 
coordinate system is located at the center of the configuration as 
indicated by the limits of integration in the following expression for 
the double moments (hereinafter referred 1o as moments): 

1 b 

Mp k = | f 2 l f 2 h t ( x .y) xk y Pd y (9) 

2 2 


where t(x,y) is the thickness distributicn of the given configuration 
and b/2 is the wing semispan. In equation (8), gnpk are constant 
coefficients given by 



^npk ^nkp l 


(n-p-k- 2 ) [ (l/2) (n+p-ik) ] I 2 k 2 P 
[ (l/2) (n-p-t - 2 ) ] 1 kl pi 


for even values 
of (n-p-k) 


|zero otherwise 


(10) 


Upon substitution of equation ( 7 ) in equal ion ( 3 ) and use of the result 
in equation (2), the integration with respect to 6 can be accomplished 
in closed form. In this manner the zero-] if t wave drag at any supersonic 
Mach number is expressed in a series in powers of p with coefficients 
which are functions of the moments; that is, 

00 

D(P) = q nDjr ( p ) 

n =2 


(ID 
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where 


Dn(P) 


1 

2it 



Un(M)] 2 as = 4 




J ^ Lnp(P c ° s 


l_P=o 



( 12 ) 

It is practical to determine the moments and make subsequent calcu- 
lations of the zero-lift wave drag with ordinary desk calculators. Con- 
sequently, calculations by the moment-of-area rule are relatively simple 
compared to those by the supersonic area rule. However, as a result of 
using a constant length for each equivalent body of revolution, the 
moment -of -area-rule series expression for the wave drag does not converge 
as rapidly as does that for the supersonic area rule. Therefore, in 
order to obtain similar solutions by both methods more terms of the 
moment-of-area-rule drag equation must be employed. As pointed out in 
reference 5 for the moment-of-area rule and demonstrated in reference k 
for the supersonic area rule, the convergence of the series can be 
improved for each method by calculating the total zero -lift wave drag 
as the s urn of the wave drag of each component alone plus the interference 
between components. For this procedure the smallest allowable length of 
each component or pair of components is employed in calculating the 
separate wave-drag contributions. 


Development of the Method 


The present method is developed by combining the concepts of the 
moment-of-area rule at Mach number 1 and the supersonic area rule. This 
combination consists merely of calculating the zero-lift wave drag by 
the moment-of-area rule at Mach number 1 for each equivalent body of 
revolution. For Mach number 1, equation (12) depends only upon the 
moments for p = 0, and can be written 


where 


D n (M=l) 




k=o 


(13) 


(14) 


Hence, for Mach number 1, equation (ll) becomes 
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D(M=1) = i q ) nDn(M=l) (15) 

n=2 


In view of equation (13), the zero-lift wave drag, at Mach number 1, of 
the equivalent bodies of revolution required by the supersonic area rule 
can be expressed as a function of the moments and lengths of the equiva- 
lent bodies as 


where 


D(M) = 




U(M)J 


nDn(M) 


n-2 


(16) 


and 


Dn(3j0) - [Lqo(Pj ( 0] 


n-2 


^ &nok 


-*k 


k=o 




MokO,0) 


(17) 


( 18 ) 


In equation (18) , the quantities Mq^(P, 0) represent the moments of the 
equivalent bodies of revolution calculated about only the spanwise axis 
through the midpoint of the length of each equivalent body (see eq. (9) 

for p = 0). From a comparison of equations (3) and (l6) it is clear 
i4 

D n (P,0) I s merely the square of the Fourier coefficients . 

Even though the concepts of the moment -of -area rule are employed to 
evaluate the Fourier coefficients, it should be noted that the actual 
length of each equivalent body, z(P,9), is considered (see eqs. (l 6) 
and (l8)). Finally the zero-lift wave drag of a given configuration is 
obtained by use of the results of equation (l6) in equation ( 2 ). 

Convergence problem .- Use of the moments to evaluate equation (17) 
permits zero-lift wave-drag calculations to be made with ordinary desk 
calculators. Experience has demonstrated, however, that the magnitude 
of the moments are such that six significant figures are required to 
obtain accurate results. As a consequence, it Is not practical to make 
calculations for more than nine terms of the series (N = 9) using desk 
calculators. In general, as few as nine terms will not provide adequate 
convergence of the series (see ref. 2 ). However, as suggested in ref- 
erence 5 and demonstrated in reference 4, for any given number of terms 
the convergence of the series can be improved for multiple component 


that k 


Lz(p,0). 
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configurations by calculating the total zero-lift wave drag as the sum 
of that for each component alone plus the interference between components. 
As a result of the ease in finding Mok(P,0) for the various components 
of a given configuration compared to finding the area, the first deriva- 
tive, or the second derivative of the area distribution of the components, 
the present method is more easily adapted to total zero-lift wave-drag 
calculations by this technique than are the methods of references 2, 3 , 
and 4. 

Dimensionless drag equations .- Zero-lift wave-drag calculations by 
the present method are most conveniently made using a dimensionless 
coordinate system. Consequently, the dimensionless quantities and zero- 
lift wave-drag equations for the present method are presented in 
appendix A. In appendix A the relationship between the dimensional and 
dimensionless zero-lift wave-drag equation for airplane-type configura- 
tions is shown to be 


D(P) 



(19) 


where 


5 



1 

2jc 



D 




c 


-.4 


Lz(pb/c, e)J 


CO 



n=2 


nDn 



Dq 





■^no 



n-2 


^nok 


Jc 


k=o 


U(pb/c, e)J 


^ok 



( 20 ) 


( 21 ) 


( 22 ) 

(23) 


The subscript 1 in equation ( 19 ) indicates properties of the exposed 
wing measured at or relative to the average chord at the wing-body 
juncture . 


Basic data, moments and lengths . - Success of the present method 
obviously depends upon a convenient technique for finding the moments 
and corresponding lengths of the various components of the equivalent 
bodies of revolution of a given configuration. The area distributions 
of body-of-revolution components of a configuration are independent of 



Mach plane orientation. Consequently, the moments of hody-of -re volution 
components are identical for all equivalent bodies of revolution. Further- 
more, within the slenderness requirements of the basic methods, the normal 
cross-sectional area distribution of body-of -re volution components, 

S(x), and the actual length of the body, Z, can be employed. Hence 
for bodies of revolution 


Mok = % S(x)^dx 

2 

or in dimensionless coordinates (see appendix A) 

Mok=§ J S(m*AS 


(24a) 


(24b) 


The area distributions and lengths of wi ng and tail components, on 
the other hand, change with Mach plane orientation. The concept of 
"sheared" configurations, which is briefly reviewed in appendix B, can 
be employed to find the area distributions of these equivalent -body 
components. In reference 4, the contribution of wing and tail components 
to the area distribution of the equivalent oodles of revolution was 
determined as a function of the tangent of the sweep angle of the 
50-percent chord line of each sheared half plan form. For a given Mach 
number and Mach plane orientation, this single parameter (K, in dimen- 
sionless coordinates) relates each sheared naif plan form to the proper 
equivalent body by the expressions (see appendix B) 


K U = K o “ P c cos 3 1 
K L = K o + P | cos 3 


(2‘j) 


where K c is the tangent of the sweep angle of the 50 -percent chord line 
of the given dimensionless plan form; Kp and K L represent the tangent 
of the sweep angles of the sheared dimensionless half plan forms in the 
positive and negative T) directions, respectively. As indicated by the 
subscripts U and L, these half plan forms nereinafter will be referred 
to as the upper and lower sheared half plan forms, respectively. The 
single parameter K is also used in appendix C of this report to deter- 
mine the moments and lengths of the sheared half plan forms. From 
equation ( 25 ) it can be seen that the quantities in equations ( 20 ) to ( 23 ) 
which are a function of Pb/c and 9 are a function of (pb/c)cos 0. 
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Procedure for Applying the Method 


General configurations . - The total zero-lift wave drag of any general 
configuration is calculated by the present method as the sum of the wave 
drag of each component alone plus the interference between all pairs of 
components, as was described in reference 4. For instance, for a config- 
uration with two components, E and H, the total zero-lift wave drag can 
be written 

5 t (p I) = Dg (p F) + Dg (p |) + Deh (p 1) (26) 


where the first two terms represent the zero-lift wave drag of each 
component alone and are evaluated by equations (20) to ( 23 ) using the 
subscripts E and H in each of the equations. The last term represents 
the mutual interference between the two components and is evaluated by 
equation (20) using the subscripts EH and the following equations: 


I P 0 ) = 2 


n 4 oo 


L Z EH( (3 " b / c ^ 0 )^ 


I 

n=2 


nD EH n ( P 9 


(27) 


' D EB n d b 9 ) = A no (p d ) A n n ( P 9 


( 28 ) 


b 


V P c' 0 


n-2 


^nok 




k=o 


L^EH ( P"^/ G > 0 )J 


^ok c> 0 


(29) 


n -2 


%no l P c i 0 


&nok 


k=o 


L^eh^^/ 0 * 0 )J 


e s 


') 


(30) 


It is emphasized that all moments are calculated about the midpoint of the 
smallest allowable length of each component or pair of components for which 
the drag is to be computed. That is, in calculating the zero-lift wave 
drag of each component alone the smallest allowable length is obviously 
the separate length of each component. For the interference term, however, 
ME Q k(P^/ c > 0 ) % 0 ^(P t) / c ^ 0 ) are calculated about the midpoint of the 

total combined length of any two components . 
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Complete wings.- Assume that the total zero-lift wave drag is to 
he calculated at "a "specified supersonic Mach number, M ma x, for a complete 
wing with 0° sweep of the pO -percent chord line (K 0 = 0) as shown by the 
solid lines in sketch (b). With regard to the airfoil section it should 

be recalled that at super- 
sonic Mach numbers for which 
the Mach waves of the flow are 
inclined ahead of the wing 
leading edge the theory is 




applicable for arbitrary air- 
foil sections; that is, both 
round- and sharp-nosed 
sections. On the other hand, 
at higher supersonic Mach 
numbers for which the Mach 
waves of the flow are 
inclined behind the wing 
leading edge, the theory is 
applicable to wings with 
KL mnu ! only a sharp-nosed airfoil 


section. For the problem 


assumed above consider a 


Sketch (b) 


wing with sharp edges. 


The calculations are initiated by a determination of the limiting 
values of Kp . and from equation (25) for 0=0. The corre- 

sponding sheared configuration is shown in sketch (b) by the dashed lines. 
(Note that |Kp . n |and | are greater than 1. If a round-nosed 

airfoiL section had been considered, 
these Limiting values would be 
restricted to values less than 1.) 

The quantity Dp(pb/c, 0) or Dy(K) 
is determined for one wing panel 
alone at arbitrary but sufficient 
values of K to define the curve 
shown in sketch (c). Addition of the 
two expressions in equation (25) 
relates the upper and lower wing 
panels in ordered pairs 

% + K L = 2 K 0 (31) 

^ use e lotion (31) > the data 
for o *8 s-f --+■ k l for ose <-f-— «-| in sketch (c) are used to make 

.. individual or combined plots of the 

b v r -7 H zero-lift wave drag of each wing 

k 0 for O40s -J-l Z Kl for 0 s8s -f panel alone as a function 

Sketch (c) 


D W (K) 
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of (3b/c)cos 0 (see sketch (d)). Since 
sheared configurations for complete 
wings are symmetrical in jt/2, positive 
values of (pb/c)cos 0 are more con- 
veniently used for the abscissa for 
ft/2>0>O; hence, in sketch (d) 

0 <(Pb/c)cos 0 < KL max . The inter- 
ference between wing panels is calcu- 
lated as a function of Kl but plotted 
directly as a function of (pb/c)cos 0 
as shown in sketch (d). The total zero- 
lift wave drag or the separate contri- 
butions can be obtained from the proper 
curve in sketch (d) by replotting the 
results against 0 for jt/2 >0 >0. 

The resulting plot is then integrated 
graphically and this result is divided 
by jc/ 2 to obtain a dimensionless value 
of the zero-lift wave drag at M max or 
(pb/c) max . The drag in pounds is, of 
course, given by equation (19)* It 
will be noted that for a complete wing 
(b^yb) and (t^/t) in equation ( 19 ) are 
unity. 



/3 -4- COS Q 


Once the zero-lift wave drag is 
calculated at M max one has inher- 
ently obtained the data from which the drag can be calculated for all 
1 < M < Mmax. For example, for 1 < M 2 < M max equation (25) yields a new 
maximum value of (pb/c)cos 0 as Kq JO (see sketch (d)). Hence, for 
(pb/c) 2 = Kl 2 , the portion of the curve in sketch (d) for 
0 <(pb/c)cos 0 < Kl 2 can replotted against 0 for jt/2 > 0 > 0, 
graphically averaged, and divided 
by jt/2. In this manner the zero- 
lift wave drag can be obtained for 
all 1 < M < Mmax or 
o <(pb/c) <(pb/c )max- 

The basic data for one com- 
plete wing panel shown in 
sketch (c) can be used to cal- 
culate the drag of each wing 
panel of the complete wing shown 
in sketch (e). This wing 
differs from that of sketch (b) 
only by the sweep angle of the 
50-percent chord line, K 0 = 1. 



Sketch (e) 
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6(0|,e) 


From use of K 0 = 1 in equation ( 31 ) it is a simple matter to obtain 
D(pb/c, 0 ) from the data of sketch (c) in new ordered pairs for 
K c < Kp < Kp max • With the aid of equatior. (25) these results for the 
wing panels alone are plotted for 0 <(pb/c)cos 0 <(KL ma x - 1 ) as shown 

in sketch (f). Because the separate wing 
panels of the sheared configurations of the 
K 0 = 1 plan form combine in different ordered 
pairs from those for the K 0 = 0 plan form, 
the interference between wing panels will be 
different for each plan form. Consequently, 
additional calculations must be made for the 
interference between wing panels for the 
K 0 = 1 plan form. It should be noted that 
the range of Mach numbers covered by the 
value ([3b/c)cos 0 = (Kp ^,^ - K 0 ) will be 
different for each plan form. 

Airplane-type configurations . - For a wing- 
body combination the zero-lift wave drag of the 
body alone is obtained from one set of calcula- 
tions, since, as noted earlier, S(x,|3,0) for 
each equivalent body of this component is 
assumed equal to S(x) . The wing-body inter- 
0-i-cose ference is calculated as a function of K or 
(pb/c)cos 0 and averaged graphically for 
selected values of pb/c. The exposed wing is 
treated in the same manner just described for 
a complete wing. In fact, in some instances, 
the basic data, Dw(K), for the wave drag of one panel of a complete wing 
can also be used to calculate the wave dra^ of each exposed wing panel 
alone of similar wings mounted on a body. This is possible whenever the 
dimensionless plan form of the exposed win 3 panel and of the complete wing 
panel are identical. As noted in appendix A, exposed wing panels are made 
dimensionless with respect to conditions at the wing-body juncture. Hence, 
D(K) for a complete wing panel can be used to calculate the wave drag of 



Sketch (f) 


an exposed wing panel alone for all straigit-line plan forms with zero 
taper ratio. For example, the wave drag o? each of the exposed wing panels 
of the wings of sketches (b) and (e) mount jd on a body can be computed as 
previously described using the data of sketch (c) for complete wings. 

The interference between exposed wing pane us must be obtained by addi- 
tional calculations since, for a given K, the combined length of both 
sheared wing panels is different for the e cposed and complete dimensionless 
wings. In computing the drag in pounds the proper values of (bi/b) 
and ( t x /t ) must be inserted In equation (li»). 


Rather than present zero-lift wave drag as a function of Mach number, 
it has been found convenient to plot reduced drag as a function of pb/c. 
One can express equation (19) in coefficient form as 
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The wing aspect ratio A is given by b 2 /S and the airfoil-section 
maximum thickness ratio, r, is given by t Q /c. Hence equation (32) can 
be written 



Since b/c is proportional to aspect ratio, equation ( 33 ) I s merely an 
expression of the linearized supersonic -flow similarity parameters; that 
is, the reduced wave-drag coefficient, Cd/At 2 , is a function of only PA. 
Thus, by making the calculations in dimensionless coordinates for a single 
configuration, one inherently computes the zero-lift wave drag for an 
entire family of related configurations . 

Systems of bodies of revolution .- The principles of the foregoing 
analyses of airplane-type configurations can also be applied to systems 
of bodies of revolution. One first renders the system dimensionless as 
outlined in appendix A. The zero-lift wave drag can be calculated as 
noted earlier for each body alone. The interference between all pairs 
of bodies can be obtained in a manner similar to that described for 
complete wing panels merely by using the parameter 6 (see appendix B) 
rather than K. A technique for finding the equivalent-body length and 
the moments of any two sheared bodies about the midpoint of their combined 
length as a function of b can be obtained by following the same 
technique outlined in appendix C for the interference between wing panels. 
It should be noted, however, that the sheared configurations of a pair of 
bodies of revolution are symmetrical in re rather than rc/2 as in the 
case of wings and horizontal tails. If an electronic computing machine 
is available, more accurate results for systems of bodies of revolution 
can be obtained more rapidly by the method of reference 4 than by the 
present method. 


EVALUATION OF THE METHOD 


The accuracy and computational time required of the present method 
will be evaluated by comparing zero -lift wave -drag solutions for several 
simplified configurations computed by the present method with solutions 
obtained by other methods. Sample calculations for the solutions obtained 
by the present method are presented in appendix D. 
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Accuracy 


Single body of revolution .- Consider the very special case of a Sears- 
Haaek hody of revolution at M = 1. As shown by slender-body theory in 
reference 6, this body has minimum zero-lift wave drag at sonic speed 
for a given length and volume. The area distribution of this body is 
given by 

3 

S(x) = S 0 1 - (h^j 2 ( 3 k) 

where S Q is the maximum frontal area of the body and l is the total 
body length. 

As shown in appendix A of this report the zero-lift wave drag of a 
system of bodies of revolution can be expressed in terms of the dimension- 
less drag as 



According to reference 6 the slender-body-theory value of D(|3d/z) at 
M = 1 is 9/8. Zero-lift wave-drag calculations by the present method yield 
exactly 9/8 (see appendix D) . It is interesting to note that the series 
expression for the zero-lift wave drag converges absolutely for the term 
n = 1. The contribution of all other terns is identically zero. This 
result further demonstrates that a Sears -Haack body has minimum zero -lift 
wave drag at M = 1 for a given length anc volume. That is, M QQ is 
proportional to the volume of a configuration (see eq. (9) for p = 0) 
and this moment alone appears only in the n = 2 term (see appendix A). 

Pair of bodies of revolution. - 
Dimensionless zero -lift wave -drag 
values have been calculated at M = 1 
for the pair of Sears -Haack bodies 
shown in sketch (g). Values computed 
by the present method (see appendix D) 
are tabulated below together with 
analytical values obtained by slender- 
body theory from reference 4. The 
total zero-lift wave drag calculated 
by the present method is accurate to 
within 9.3 percent of the slender-body- 
theory solution of reference 4. As 
noted earlier for a single body of 
revolution, values for the bodies 



Sketch (g) 
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BL f3(d/z ) ] for 

bodies in sketch 

(g), M = 1 

Drag component 

Present method 

Reference 4 

Body E alone 

1.125 

1.125 

Body H alone 

1.125 

1.125 

Interference 

.204 

.340 

Total 

2.454 

2.590 


alone by the present method are in exact agreement with the analytical 
solutions of reference 4 but the interference values agree to within only 
40 percent. This large discrepancy for the interference calculation 
results from the fact that the moments for this calculation must be 
calculated about the midpoint of the combined length of both bodies , a 
length longer than that of the individual bodies. As noted in refer- 
ence 5; when this occurs the contributions to the drag of the terms 
for n > 9 (higher harmonics) become significant. 

Complete wings . - Zero-lift wave-drag calculations have been made for 
two families of complete wings with biconvex airfoil sections and a 
constant thickness ratio from root to tip. One family has a diamond 
plan form (K 0 = 0) and the other a triangular plan form (K 0 = 1) as 
shown in sketches (b) and (e), respectively. The linearized theory 
zero-lift wave-drag solutions for these wings have been obtained in 
reference 7 hy a stepwise integration of the pressures over the entire 
wings. The zero-lift wave-drag results from reference 7 and those 
obtained from calculations by the present method for N = 9 (see appen- 
dix D) are presented in figure 1 in reduced drag-coefficient form. In 
order to evaluate the present method relative to the method of the super- 
sonic area rule and the moment-of-area rule, zero-lift wave-drag solutions 
for both families of complete wings were calculated by both latter methods, 
the two methods upon which the present method is based. These latter 
results are also presented in figure 1. Results for the supersonic area 
rule are usually obtained by the method of reference 2. In this report, 
however, the results for N = 9 were obtained by the technique of the 
present method by calculating the total moments of each equivalent body 
about the midpoint of the total length of each sheared configuration. 

By this approach the equivalent-body zero-lift wave-drag equations for 
the present and supersonic -area-rule methods can be shown to be identical 
theoretically. Numerical calculations by both methods (N = 9) for several 
equivalent bodies agreed to within one half of 1 percent. Results for 
the moment-of-area rule were calculated for N = 9 by the alternative 
drag formula given in appendix C of reference 5. As indicated in figure 1, 
for the range of values of (Bb/c for which data are available, the present 
method yields results which differ from the linearized- theory solutions 
of reference 7 between 1 and 6 percent. The data of figure 1 also 
indicate that the solutions obtained by the present method are more 
accurate than those obtained by the other two methods. Furthermore, 
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solutions obtained by the supersonic area rule are more accurate than 
those obtained by the moment -of -area rule. The more accurate solutions 
by the present method are a result of calculating the drag of the wings 
by parts. Since the solutions obtained by the three "series” methods 
were calculated for N = 9 i n each case, the relative accuracy of the 
solutions obtained by these three methods is associated with the relative 
rapidity of the convergence of the series expression of the zero-lift 
wave drag for each method. The relative convergence of the series for 
the three methods will be discussed subsequently. 

Airplane-type configurations .- The procedures for applying the present 
method to airplane -type configurations have been discussed earlier in 
this report and in appendix A. Although no calculations have been made 
to evaluate the accuracy of the present method for calculating the zero- 
lift wave drag of such configurations relative to other methods, it is 
felt that differences in such solutions would be comparable to those 
discussed above for complete wings. However, it should be noted that, 
as discussed in appendix A, some errors, not present in the case of 
complete wings, might be introduced into the calculations by the present 
method as a result of assuming a constant wing station for the wing-body 
juncture when computing the moments of the exposed wing panels. These 
errors will be more significant in cases of severely indented bodies and 
in cases where the body radius is of the same order of magnitude as the 
wing thickness at the wing-body juncture. 

Convergence problem .- As mentioned in reference 5* in the Fourier 
series analysis of the slope of the area distribution of a configuration 
(a type of analysis basic to the present method, the supersonic area rule 
and the moment -of -area rule), the series coiverges most rapidly when the 
smallest allowable length of the configuration is used in the wave -drag 
analysis. The total zero-lift wave-drag results shown in figure 1 for 
the moment-of-area rule used a constant total length as specified in 
appendix C of reference 5- The supersonic area rule employed the total 
length of each individual sheared complete *ang. Use of the individual 
total lengths with the supersonic area rule results in a marked increase 
in the rapidity with which the series converges. This is indicated in 
figure 1 by the larger values of the wave drag obtained by the supersonic 
area rule as compared to those obtained by bhe moment-of-area rule. The 
total zero-lift wave drag was calculated by the present method as the 
sum of the wave drag of each wing panel aloie plus the interference between 
wing panels. The individual length of each wing panel of each sheared 
configuration was used for the drag calculations of each wing panel alone 
and the combined length of both sheared wings panels was used for the 
interference calculations. In this manner the present method provided 
an additional increase in the rate at which the series converges. Thus, 
the values of the wave drag in figure 1 for the present method are greater 
than those obtained by the supersonic area rule and the moment-of-area rule. 
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Regardless of whether the total zero-lift wave-drag calculations are 
made using the smallest allowable length or the total length of each 
sheared configuration, the interference calculations in both cases are 
made using the total length. On the other hand, calculations of the 
wave drag of each wing panel alone are made using either the length of 
each sheared wing panel or the combined length of both sheared wing panels. 
Hence, the magnitude of the increased accuracy of the total zero-lift 
wave -drag calculations afforded by the use of the smallest allowable 
length depends upon the ratio of the length of each sheared wing panel 
alone to the combined length of both sheared wing panels. For convenience 
this ratio is designated by e. The smaller the value of e compared 
to unity the greater wall be the improvement in the drag calculations as 
a result of using the smallest allowable rather than the total length. 

This fact is demonstrated by the data in figure 1. Up to a Mach number 
for a sonic leading edge for the diamond plan form ((3b/c = 1 in 
fig. 1(a)) e = 1 for both wing panels of all sheared configurations; 
hence, the single curve for the results of the present and supersonic- 
area-rule methods in this region. For P'b/c > 1, e is less than 1 
for only some of the sheared configurations and is also the same value 
for both 'ring panels. In this region, therefore, only a very slight 
improvement in the accuracy of the solutions obtained by the present 
method is realized compared to solutions obtained by the supersonic area 
rule. For the triangular plan form (fig. 1(b)) e = 1 only for Q = jt/ 2 
or for pb/c = 0. For all other sheared configurations e < 1 and has 
different values for each wing panel. Thus, the accuracy of the solutions 
obtained by the present method for the triangular plan form is substan- 
tially improved compared to the accuracy of the solutions obtained by the 
supersonic area rule. 

In view of the benefits demonstrated above resulting from the calcu- 
lation of total zero-lift wave drag by parts, the question naturally arises 
as to the possibility of using this procedure in conjunction with the 
supersonic area rule. The method of reference 2 can, in fact, be adapted 
to such a procedure. This, however, is not recommended unless one is 
prepared to accept a two- to threefold increase in work load. 


Computing Time Required 


The method of reference 7 uses the pressures over the entire ring 
plan form to calculate total zero-lift wave drag. These calculations are 
necessarily tedious and lengthy. The actual computational time is not 
given in reference 7- However, from a comparative examination of the 
wave -drag equations this method obviously requires more computational 
time than the present method or the two methods upon which the present 
method is based. Furthermore, the latter three methods each employ 
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some geometric property of the configuration to evaluate the Fourier 
coefficients of a series expression for the wave drag. For these reasons 
the computational time required of present method is evaluated relative 
to only the supersonic area and moment -of -are a rules. The computational 
time in each case is based on a series for N = 9* 

The basic differences in the computational time required of the 
moment-of-area rule, the present method, and the supersonic area rule 
are a function of the time required by eacl method to obtain the respective 
geometric properties used to evaluate the Fourier coefficients. The actual 
time, of course, depends upon the complexity of each configuration. The 
total zero -lift wave -drag calculations for the diamond-plan -form wing by 
the moment-of-area rule required 20 man -hoi rs with desk calculators. By 
this method the for only the original complete wing were required. 

No graphical averaging was required. The present method, applied with 
either the smallest allowable or the total length, required 100 man-hours 
with desk calculators. By this method the Mok(ph/c, 0) for 20 sheared 
configurations were used. The wave-drag results for these configurations 
we re graphically averaged. No calculations were actually made by the 
supersonic area rule; however, a computing time of 500 man-hours has been 
estimated for this method by the authors of reference 2. This large 
estimate was based primarily upon the determination, with desk calculators, 
of the area distribution of the 20 sheared configurations and secondarily 
upon a graphical averaging of the wave -drag results. 

The basic data for the wave drag of each sheared wing panel alone 
of the diamond plan form was used in the manner described in the METHOD 
section to obtain the wave drag of each wii.g panel alone of the triangular 
plan form. The only additional calculatiox. required to obtain the total 
zero-lift wave drag of the triangular plan form by the present method was 
that for the interference between wing panels and, of course, the graphical 
averaging of all results. In this manner only 30 man-hours were required 
by the present method to obtain the results for the triangular-plan-form 
wing as compared to 100 man-hours for the clamond-plan-form wing. Similar 
to the case of the diamond-plan-form wing the moment-of -area-rule calcu- 
lations required 20 man-hours for the triangular plan form. Since the 
method of reference 2 does not normally yield the wave drag by parts, 
zero-lift wave-drag calculations by this method would also require an 
estimated 500 man-hours for the triangular-plan-form wing. 

The 500 man-hour computation for the supersonic area rule was made 
without a knowledge of reference 4 or the present report which both utilize 
the parameter K to organize the basic data for wings. It is reasonable 
to expect, therefore, that use of the parameter K, and/or mailing calcu- 
lations for related wing by parts, would result in a lower estimate of the 
computing time required of the supersonic area rule. Furthermore, the 
differences in computing time, of course, would be substantially reduced 
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if the complete procedures for all methods were programed on electronic 
computing machines. It is emphasized, however, that the time required 
to program each method would be different and is unknown at present. 


CONCLUDING REMARKS 


A method has been developed for calculating zero-lift wave drag 
by combining the concepts of the supersonic area rule and the moment- 
of-area rule. The accuracy of the method has been evaluated by comparing 
total zero-lift wave-drag solutions obtained for several simplified con- 
figurations by the present method with solutions given by slender-body 
and linearized theory. The accuracy and computing time required of the 
method have also been evaluated relative to the two methods upon which 
the present method is based. The following remarks are warranted as a 
result of the development and evaluation of the present method. 

Results are obtained by the present method from a series expression 
for the zero-lift wave drag which can be evaluated with an ordinary desk 
calculator. The total zero-lift wave drag of a configuration is calcu- 
lated as the sum of the wave drag of each component alone plus the inter- 
ference between components. Total wave-drag results for several simplified 
configurations demonstrated that, by calculating the total drag by parts, 
nine terms provide adequate convergence of the series. Total zero-lift 
wave -drag results obtained by the present method differed from solutions 
given by slender-body and linearized theory by less than 6 percent for 
all configurations investigated. For nine terms of the wave -drag series, 
these results represent an accuracy greater than that obtained by the 
moment -of -area rule and greater than or equal to in some isolated cases, 
the results obtained by the supersonic area rule. The greater accuracy 
of the present method is a direct result of calculating the total zero- 
lift wave drag by parts. 

The computational time required by the present method was five times 
as great as that required by the moment -of -area rule and one-fifth of 
that required of the supersonic area rule.. The shorter computing time 
of the moment-of-area rule, however, is of no real consequence in view 
of the poor accuracy of results obtained by this method. 'The differences 
in computing time, of course, would be substantially reduced if the 
complete procedures for all methods were programed on electronic computing 
machines . 

As a result of making the calculations in dimensionless coordinates 
the basic data obtained for a configuration at a specified supersonic 
Mach number can be manipulated to obtain the total zero-lift wave drag 



of the given configuration for all lower supersonic Mach numbers with a 
minimum of additional calculations. Furthermore, the basic data for a 
given configuration which is used to calculate the wave drag of each wing 
panel alone can he manipulated to provide the wave drag of each wing panel 
alone of an entire family of similar wings 'which differ only in sweep. 


Ames Research Center 

national Aeronautics and Space Administration 
Moffett Field, Calif., Jan. 19 , 1959 



APPENDIX A 


DIMENSIONLESS ZERO -LIFT WAVE -DRAG EQUATIONS 


In order to take advantage of the decrease in the number of parameter 
resulting from similarity considerations, and to facilitate the calcula- 
tions, the quantities defined in the METHOD section can he made 
dimensionless . 


AIRPLANE -TYPE CONFIGURATIONS 
Complete Wings Alone 


As shown in appendix B of reference p, all x coordinates can be 
divided by the half -chord of the wing at the vertical plane of symmetry, 
c/2; all y coordinates can be divided by the semi span of the complete 
wing, b/2; and all z or thickness coordinates, t, can be divided by 
the wing maximum thickness at the vertical plane of symmetry, t Q . In 
this manner the dimensionless quantities and their relationship to the 
corresponding dimensional quantities can be defined as follows: The 

dimensionless Cartesian coordinates are defined by 
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b c 


2y 
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r- or t(|,t)) 


t( x >y) UQti) 


(Al) 


■The moments for the equivalent bodies of revolution can be written with 
the aid of equation (9), for p = 0, as 


b Z(p,9) 
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With the definitions of equation (Al) , equation (Al) is made dimensionless 
as follows : 
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Finally, the relationship between the dimensional and dimensionless zero- 
lift wave drag for an equivalent body of revolution of a complete wing is 
given by (see eq. (l6)) 
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Inspection of equation (2) 


indicates that for a given configuration 


D(0) 






(A8) 


where 
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Wings in the Presence of Bodies 


For a given configuration the present method requires separate calcu- 
lations of the wave drag of each component alone. Hence, the r\ limits 
of integration in equation (A3) for one panel of an exposed wing (or tail) 
w T ould "be from the wing-body juncture, rj = 1^(5), to the wing tip, r\ ~ 1. 
Calculation ot the moments is obviously simplified if the rj limits of 
integration are constant from 0 to 1. 'This is accomplished by assuming 
a constant wing station for the wing-body juncture which is an average 
wing station of the actual juncture and by rendering the resulting exposed 
wing panel dimensionless with respect to the conditions at this averaged 
wing-body juncture. Use of the averaged wing-body juncture rather than 
the contour of the actual wing-body juncture has been found to have a 
negligible effect upon the drag calculations except in cases of severe 
body indentations and in cases where the body radius is of the same order 
of magnitude as the wing thickness at the wing-body juncture. Thus 
the x coordinates can be divided by the half-chord of the wing at the 
juncture, c 1 /2, the y coordinates can be divided by the semispan of the 
exposed wing, b 1( /2, and the thickness coordinates can be divided by the 
wing maximum thickness at the juncture, t^_. In this manner equation (A8) 
can be written 



Most practical wing (or tail) plan forms have straight leading and trailing 
edges; hence, b 1 /c 1 = b/c. It is convenient, therefore, to write equa- 
tion (A10) in a form similar to that of equation (A8); that is, 
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SYSTEMS OF BODIES OF RESOLUTION 
Single Body of Re volition 


The moments of a single body of revolution are determined from the 
irca distribution of the body as given by equation (24a); that is, 
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It is convenient, therefore, to render a body of revolution dimensionless 
by dividing the x coordinates by the half -body length, lj 2, and by 
dividing the area distribution by the maximum frontal area of the body, S Q « 
Thus 
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In this manner the dimensionless moments of e .matiort (Al-;) can be defined 
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Equations for I^q, iL, and D for a single body of revolution are obtained 
from equations (A4) , (At), and (A7)j, respectively, by replacing t 0 (b/h) 
and c/2 by S G and lj 2, respectively. A single body, of course, is 
independent of |i and 9 . Hence, from equation (A6) or (Ab) the relation- 
ship between the dimensional and dimensionless zero-lift wave drag for a 
single body of revolution is given by 
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Two or More Bodies oT Revolution 


Individual dimensionless coordinate systems are used to calculate 
the wave drag ox ea.cn tony alone . To calculate the interference between 
pairs 01 oodles ox revolution, however, it has been found more convenient 
to render each pair of bodes dimensionless with respect to the length 
and maximum frontal area of only one of the bodies. Usually the largest 
oody is selected. The dimensionless quantities L no , Dn, and D for 
interference calculations between pairs of sheared bodies (see appendix 13) 
are obtained in the manner described above for a single body. In addition, 
quantities which are a function of (pb/c, 0) are replaced by (pd/z , 0), 
where d is the lateral spacing of the bodies, as shown for example in 
sketch (g), and l is the length of that body upon which the dimension- 
less quantities arc cased. 


For the convenience of those who may use the present method the 
Lno(P° d fl) equation (A4) are listed below for n up to and includ- 
ing ta 'Fh.e parameters pb/c and 0 have been omitted from the notation 
in the interest of simplicity. 
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APPENDIX B 

CONCEPT OF SHEARED CONFIGURATIONS 


With the supersonic area rule, one calculates the zero-lift wave 
drag of a configuration as the average of that of a series of equivalent 
bodies of revolution. The normal cross-sectional area distribution of 
these equivalent bodies is obtained as the frontal projection of the area 
distribution intercepted on the given configuration by a set of parallel 
oblique planes tangent to the Mach cones. In reference 8 a new concept 
was first introduced which permits the use of the normal cross-sectional 
area distributions of a series of sheared configurations which, within 
the slenderness requirements of the supersonic area rule, are considered 
equal to those of the equivalent bodies obtained in the manner just 
described. In reference 4 the concept of sheared configurations was 
employed as a simplified technique for fincing the contribution of wings, 
tails, and pairs of bodies of revolution to the area distribution of the 
equivalent bodies of revolution. The present method employs the concept 
of sheared configurations as a simplified technique for finding the 
moments and lengths of the various components of the equivalent bodies 
of revolution. Therefore, a brief review of this concept is presented 
below for wings, tails, and pairs of bodies of revolution. 

WINGS AND TAI U 
Wings 


Complete wings .- A complete wing is shown in dimensionless coordinates 
by the solid lines in figure 2(a). This wing is made dimensionless with 
respect to conditions at the vertical plain of symmetry (see appendix A). 
The traces in the plane (see fig. 2(a)) made by the intersection 

of the oblique planes and the £,tj plane are inclined to the r\ axis 
at an angle ^ so that >jr = tan" 1 [ (pb/c ) cc s 0]. Hence, the traces 
through the 50-percent chord at the wing tip and the £ axis (see 
lines XX and YY in fig. 2(a)) define a lcngitudinal distance 
| (pb/c)cos 0|. In accordance with reference 4 the sheared configurations 
are formed by shearing each element of the upper wing panel a distance 
| (f3b/c)cos 0 \r) forward and shearing each clement of the lower wing panel 
a distance | (pb/c)cos 0|q rearward for ( • < 0 < jt/2 (see the dashed 
lines in fig. 2(a)). For 0 - jt/2, the sheared and given configurations 
are identical. For it/2 < 0 < it, the upper and lower wing panels are 
sheared rearward and forward, respectively. However, in view of symmetry 
of a pair of wing panels only the values of 0 < 0 < jt/2 are considered. 
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Within the slenderness requirements of the supersonic area rule (see the 
METHOD section), the normal cross-sectional area distribution of the 
sheared configurations is identical to the frontal projection of the 
area distribution intercepted on the given configuration by the oblique 
planes. From figure 2(a) it is seen that the tangent of the sweep angle 
of the 50-percent chord line of the upper and lower sheared wing panels 
can be defined by 


K U = K o " P £ cos 0 


= K 0 + (3 cos 0 

> 


(Bl) 


Addition of the two parts of equation (Bl) yields the following relation 
for the upper and lover wing panels. 


K L + Ku = 2K 0 


(B2) 


Thus, for a given p and 0 the sheared wing panels are related to the 
proper equivalent bodies of revolution by the single parameter K. 

Exposed wings . - The exposed wing panels of a wing-body combination 
are shown in dimensionless coordinates in figure 2(b). Since these wing 
panels are made dimensionless with respect to conditions at the wing- 
toody juncture (see appendix A), it is convenient to consider a separate 
coordinate system for the upper and lower exposed wing panels . The | 
axes of these separate coordinate systems are separated in the spanwise 
direction by a distance 2r\ ± (see fig. 2(b)). The distance m is the 
ratio of the body radius at the averaged wing-body juncture to the semi- 
span of the exposed wing panels. With respect to the coordinate system 
of the complete wing, the semispan of the complete plan form is 1 + tj x . 
Hence, the trace XX in figure 2.(b) through the 50-percent chord at the 
wing tip and the | axis defines a longitudinal distance 
I (Pb 1 /c 1 )cos 0 | ( 1 + t|;l ) . The sheared configurations of the exposed wing 
panels are formed by shearing the plan form as before. As the exposed 
wing panels are sheared with respect to origin of the complete wing 
(0 in fig. 2(b)), the origins of the coordinate systems of the separate 
exposed wing panels become longitudinally separated a distance 
A = 2 | (pb 1 /c 1 )cos 0 | iqx as shown by the dashed lines in figure 2(b). 

The origin of each panel is translated, in opposite directions, a distance 
l(Pb 1 /c 1 )cos 0 | . Since b ! /c x = b/c for plan forms with straight 

leading and trailing edges, it can be determined with the aid of figure 2(b) 
that equation (Bl) also serves to identify the sheared exposed wing panels 
with the proper equivalent bodies of revolutions when the equations are 
applied relative to the individual coordinate system of each exposed 
wing panel. 
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The value of identifying the sheared wing panels with the parameter 
K rather than with 3 and 0 was first demonstrated in reference 4. In 
reference 4 it was found possible to calculate the area distribution of 
the sheared wing panels as a function of |, K, and fixed geometric 
properties of the given plan form. As shown in appendix C of this report, 
the moments and lengths of the sheared wing panels can also be calculated 
as a function of K and fixed geometric properties of the given plan form 


Tails 


'Hie shearing concept described above for exposed wings is also 
applicable to the exposed portion of tails when the tails are rendered 
dimensionless with respect to conditions at the tail-body juncture. 

The proper values of K for horizontal tails are determined from 
equation (Bl) with = 0y. For vertical tails either part of equa- 

tion (Bl) is used with 0yp = 0 W + n/2. '.lie sheared configurations for 
vertical tails are symmetrical in it rather than ir/2. 


PAIRS OF BODIES OF REVOLUTION 


A pair of bodies of revolution is shown in dimensionless coordinates 
by the solid lines in figure 3 . The configuration is rendered dimension- 
less with respect to the larger body, E, and the lateral distance between 
the bodies, d. If the origin of the coor linate system is chosen at the 
center of body E, traces in the !,t) plane made by the intersections 
of the oblique planes and the plane are inclined to the t) axis 

at an angle S so that w = tan -1 [ (pd/lg ) cos 0], Similar to the case 
of sheared wings, each element of body H is sheared a distance 
| (pd/ Zg)cos 0|. Sheared configurations for a pair of bodies of revolu- 
tion are symmetrical in it. Hence, for ) < 9 < n/2 body H is sheared 
rearward. For 9 = it/2 the sheared and given configurations are identical 
For n/2 < 0 < n body H is sheared forward as shown by the dashed lines 
in figure 3. For a given 3 and 0 the sheared configurations can be 
related to the proper equivalent bodies o? revolution by the single 
parameter & which indicates the longitudinal separation of the centers 
of the sheared bodies; that is, 


5 = L + P i c:s 0 (B3) 

where 5 0 is the longitudinal distance, Ln dimensionless coordinates, 
between the centers of the given bodies ( jee fig. 3)- 
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APPENDIX C 

MOMENTS AND LENGTHS FOR SHEARED WINGS 
AND TAILS 


Zero-lift wave-drag calculations by the present method require a 
knowledge of the moments and lengths of the various components and pairs 
of components of a series of equivalent bodies of revolution (or sheared 
configurations). As mentioned in the METHOD section, the moments and 
lengths of slender-body-of -re volution components alone are relatively 
simple to determine because the area distribution of the sheared com- 
ponents are independent of Mach plane orientation. The lengths of pairs 
of these components are also relatively simple to determine, since, as 
shown in appendix B, the effect of changes in 0 on the sheared config- 
urations of such components is merely one of translation of one body 
(see fig. 3)- The area distribution of sheared wing and tail components, 
on the other hand, is dependent upon 0. Consequently, the moments and 
lengths required are more difficult to determine. Therefore, this 
appendix presents a technique for determining the moments and lengths 
of wings and tails required by the present method to calculate the zero- 
Tift wave drag of each panel alone plus the interference between panels. 

In the following analysis, use of the term n wing(s)' 1 will be construed 
to include both wings and tails. Application of the results to tails is 
made merely by means of the proper value of K for the tails as previ- 
ously described in appendix B. 


The results of the following analyses are applicable to a large yet 
restricted group of wings. This group satisfies the following conditions: 

i. The spanwise variation of the ^O-percent chord line is linear, 

ii. The boundaries of the plan form, the spanwise variation of the 
local chord, and the spanwise variation of the thickness along lines of 
constant percent chord can be expressed analytically. 

iii. The wing has the same airfoil section at all spanwise stations; 
that is, only spanwise variations in thickness and/or thickness ratio are 
allowed. 

Condition (iii) can be relaxed to include wings with one airfoil section 
out to a discontinuity in plan form, for example, a fence or extended 
leading edge, and a different airfoil section from the discontinuity out 
to the wing tip. In such cases, each portion of the plan form is treated 
as a separate entity. 
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This appendix presents the results onLy for trapezoidal plan forms 
with taper ratios from 0 to 1. Moments ani lengths for other plan forms , 
such as an elliptic plan-form wing can be obtained by following the same 
technique presented herein and inserting tie proper analytical expressions 
for the plan form and the spanwise variation of the local chord (or 
half-chord) . 


MOMENTS AND LENGTHS FOR THE CALCULATION 
OF THE DRAG OF ONE WING PANEL 

Moments 


The sheared plan form of one exposed wing panel of a wing-body 
combination is shown in dimensionless coordinates in figure 4. As shown 
in appendix A, the moments of one sheared wing panel can be expressed as 


Mok 





Z(0b/c, 0) 


!=-Z(0b/c, 0) 



| k d| d n 


(Cl) 


It will be recalled that this definition requires that the moments be 
calculated about the midpoint of the total length of the sheared wing 
panel. It has been found expedient to calculate first the moments about 
the axis through the origin of the coordinate system shown in figure 4 
(| = 0) and then transfer these moments, designated M^^ph/c, 0 ) , to 
an axis through the midpoint of the total Length. Thus 


Lk' 




| 2 (0b/ C , 0,T)) 


|=!i(3b/c, 9,7}) 



| k d£ dr) 

(C2) 


where li(PL/c, 0,t]) and ^(Pb/c, 0,r]) are functions which define the 
leading and trailing edges of the sheared plan form, respectively. 
Equation (C2) can be simplified somewhat by replacing the combined 
variables Pb/c and Q by the single parameter K (Kp or Kp) as described 
in the METHOD section (see eq. (25)). Hence, 


M ok '(K) 



p5a(K,Ti) 

!=Si(K,ti) 


T(!,T),KH k d5 dT) 


(C3) 



As shown in appendix C of reference 4, the plan-form dimensionless 
thickness distribution for plan forms satisfying the above mentioned 
restrictions can be written 

t(£,T!,K) = 4»(Tl)T[a(5,T),K)] (C4) 

where $(tj) is the spanwise variation of the dimensionless thickness along 
lines of constant percent chord, T[a(£,q,K)] is the local airfoil section 
dimensionless thickness along lines of constant percent chord, and 
a(|,T),K) defines lines of constant percent chord. In other words, the 
function a(|,r),K) expresses any point |,t) in any sheared plan form 
(for any value of K) as the chord station of the local dimensionless 
airfoil section. The chord station is measured relative to the local 
50-percent chord as indicated by 


a(|,ri,K) = a = ilH ; -1 < a < 1 (C5) 

c(ti) 

In equation (C5) c(q) is the spanwise variation of the local dimensionless 
half-chord, which, for the plan forms considered herein, is given by 

c(ti) = 1 - Bt) (C 6) 

In equation (C 6) 

B = 1 - (C7) 

and is the taper ratio of the exposed wing panel. From equations 

(Cp) to (C7) equation (C4) can be written 

= 3>(7)T(a) (C8) 


■1 < a < 1 (C9) 

Substituting equations (C8) and (C9) and the equations for li(K,q) 
and. | 2 (K,q) (see fig. 4) in equation (C3), making a change of variable 
from ^ to a, and first performing the integration with respect to a 
yields 


whe re 


•(i,T],K) 


a .1*21 ; 
1-Br) 


Mok'(K) 



t(og) [KT)+(1-Bq)a] (l-Br))da, dr] 


a=-i 


(CIO) 



[KT]+(l-BTi)a] k = ) cJj(Kr)) k " m (l-Bri) m a n 


where 


_ k: 

°m “ *7 7 " ” 

(k-m) l ml 


lienee, equation (CIO) becomes 


Mok f (K) = > 4 IC 


K k " m § f T(a)a :n dc. r r] k “ Ifl ( 1 -B r| ) ra+ 1 ^ ( r| ) d.T| 


and since 


(1-Br]) m+1 = ) (-l) r c:™ +1 B r r 1 r 


juation (Cll) becomes 


W 

J T]— O 


c y- m Y (-dW: 


r k-m+r m 


(Cll) 


X m = § f T(a)c. m da 


(CIS) 


Mo k '(K) = ) C k K k “ m i r:1 ) (-l) r C™ + ' B 3 


T )k-m +r $(T 1 )dT} (C13) 


m=o 


r-o 
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Finally, if 


■^km 


ra+i 


(-1) V! +1 B r E 


k-m+r 


r=o 


(C15) 


the moments about the axis through the origin of the coordinate system 
of the sheared wing panel (| =0) are given by 

k 

Mok'OO = Y, (Cl6) 

m-o 


The physical interpretation of equation (Cl6) is that the moments 
are functions of the single parameter K, which identifies the sheared 
wing panel, and the geometric properties of the given wing panel. As 
indicated by equation (C 15 ), is a function of the taper ratio of 

the given wing panel and the spanwise variation of the dimensionless 
thickness of the given wing panel along lines of constant percent chord; 
that is, Ai enters through B (see eq. (C7)) and o(q) enters through 
^k-m+r (see eq. (Cl^-)). It should be recalled (see restrictive condi- 
tion (ii) above) that 0(q) is a known analytical expression; therefore, 
^k-m+r can always be evaluated with little or no difficulty. Several 
practical examples of o(q) are tabulated below for the plan forms 
discussed herein. 


Thickness 

characteristic 


Definition of constants 

Constant 
thickness ratio 

1-Bt] 

B (see eq. (C7)) 

Linear variation 
in thickness ratio 

(1-Bt}) (1-Gt)) 

n thickness ratio at rj=l 

vJT _L 

thickness ratio at rj=0 

Linear variation 
in thickness 

1-Ht] 

^ ^ thickness at q=l 

thickness at r\=0 


The quantities I m are integral functions of the dimensionless airfoil- 
section thickness distribution, t(cc) (see eq. (C12)). When T (a) is 
known analytically I m can also be evaluated with little or no 
difficulty. Generally, however, an analytical expression for t(cx) may 
not be available. In these instances the evaluation of equation (C12) 
requires special treatment. A simplified technique for finding I m for 
airfoil sections with a nonanalytic thickness distribution is presented 
in appendix E. 
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It remains to transfer the moments from the axis through | = 0 
an axis through the midpoint of the total length of the wing panel, 
the aid of sketch (h) it is seen that such an axis passes through 

7 ) 7 ? 



to 

With 


£ = A^(K) where A^(K) indicates the moment transfer distance (and 
direction relative to the axis through £ = 0) for a sheared wing panel 
alone. The analytical expression for the moments calculated about the 
£ = Av/(K) axis is 


Mok(K) =| f f TU,T!,K)[^A W (K)] k d^ dTj 

o T| ^ £ 


(C17) 


Since 


k 


U-A w (K)r = ^ (-l) S C^ k " S [A w (K)] i 
s =o 


equation (C17) can be written 
k 


M ok (K) 


(-D S cWk)] S § 


TU,T),K)| k - s d| dr) (C18) 


s=o 




By definition from equation (C3) 


fli(k-B)(K) =1 f f T(^T) ; K)^ k - S d| dT) 

J T\ J ^ 


which defines moments about the £ = 0 axis. Therefore, moments about 
the ^ = Ay(K) axis are given as a function of the moments about the 
£ = 0 axis as 



39 


k 

M 0 A K ) = X ( - l)Sc ^ A w( K )^ S i(k-s)(K) (C19) 

s=o 

Values of Ay(K) for various values of K cannot , of course, be 
determined until the total length of the sheared wing panel, 2 w ( K t is 
known. Therefore, Ay(K) is evaluated in the following paragraph in 
conjunction with the evaluation of Zy(K) . 


Lengths and Moment Transfer Distances 


The analytical expressions for the total length and moment transfer 
distance for each sheared wing panel alone have different forms depending 
upon the value of K relative to B. Sheared wing panels for ranges of 
K < -B, -B < K < B, and K > B are shown in figure 5- Consider the 
sheared wing panel for K > B shown in figure 5(c). Actually the length 
ratio [c/z(K)]y or [c/z(K)]y is of interest. From figure 5( c ) 

Zy(K) = 1 + (K+Ai) (C20) 

and 

A W (K) = (K+Ai) - | Z w (k) 
which, from equations (C7) and (C20), becomes 

A W (K) = | (K-B) (C21) 

Since c = 2, from equation (C20) 



2 

1+(K+Ai) 


which, from equations (C7) and (C21), can be written 



1 

1+A W (K) 


(C22) 



In a similar manner Ay(K) and [c/z(K)]y can be determined for the range 
of K defined by the sheared wing panels in figures p(a) and p(b). The 
results are tabulated below: 


Moment transfer distances and length 
ratios for sheared ring panels alone 

Sheared 
wing panel, 
K range 

K < - B 

-b < 1 : < b 

K > B 

AwOO 

(1/2) (K+B) 

0 

(1/2) (K-B) 

c 

1 

•L 

1 

IfObJw 

l-Ay(K) 

1+A W (K) 


MOMENTS AND LENGTHS FOR THE CALCULATION 
OF INTERFERENCE DRAG BETWEEN WING PANELS 


Moments 


Moments for the upper and lower sheared wing panels used in the 
interfere.: ice calculations are determined :.n the same general manner as 
in the case of the drag calculations for one sheared wing panel alone. 

The moments used in the drag calculations for one sheared wing panel 
alone which were calculated about the £ = 0 axis by equation (Cl6) 
car. also be used in the interference calculations to obtain the moments 
about the midpoint of the total length of the combined sheared wing 
panels. Values of Mok ? (K) are selected for the upper and lower wing 
panels in the ordered pairs indicated by equation (hi). These moments , 
in conjunction with the respective upper and lower transfer distances for 
the interference calculations, Ajy(K) and Aj-^(K), are used in equa- 
tion (C19) to obtain values of Mq^-(K) for interference. The actual 
transfer distances are, again, dependent upon the combined length of 
bo tli sheared wing panels. 


Lengths and Moment Transfer Distances 


The analytical expressions for the total length of both wing panels 
and tiie separate moment transfer distance for each ring panel used in the 
interference calculations for a sheared plan form have different forms 
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depending upon the value of K 0 relative to B and upon the value of 
Kg relative to K 0 and to 3. Sheared plan forms defining the ranges 
of critical values of K 0 , B, and Kp are shown in figure 6. 

Rather than express A(K), Aj u (K), Aj (k), and [c/z(K)] I as a 

function of Kg and Kg, it has been found convenient to make use of the 
relationship between Kg and Kg given by equation ( 31 ) and express all 
quantities as a function of Kg. For all the sheared plan forms of 
figure 6 

A(K) = tiPKl-Ku) 

but from equation ( 3 I) one can write 


A(K) = 2tu(K l -K 0 ) (C23) 

Consider the plan forms shown in figure 6 (e) for K 0 > B and 
Kg >(B+2 Kq) . From the figure 

Zj(K) = (Kg+Ai) + A(K) - (Kg- Ai) (C2k) 

Alu(K) = | lj(K) + (Ku-Ai) (C25) 

A Il (K) = (Kl+Ax) - | g(K) (Ci' 6 ) 

From equation ( 31 ) equations (Ci4) to (Cf 6 ) become 


K K) = 2(K L -K 0 +Ai) + A(K) 
Aiu(K) = i [ 2 K q +A(K) ] 

A Il (K) = | [ 2Kq-A(K) ] 


Finally, since c = 2 



1 

K L -K 0 +A 1+ | A(K) 


(C27) 

(C28) 

(C29) 


(C30) 


In a similar manner AjgK), AjgK), and [c/z(k)]j can be determined for 
the ranges of critical values of K 0 , B, and defined by the other 
sheared plan forms in figure 6 . The results are tabulated below; 
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Moment transfer distances and length ratios for the 
interference between wing panels of sheared plan forms 

Given 
plan form, 
K 0 range 

0 < K 0 < B 

0 < K 0 < 3 
and K q > -3 

0 < Kq < B 
and K 0 > B 

Sheared 
plan form, 
Kp range 

Ko < K L < B 

B* < K L <(B+.aCo) 

K l >(B+2K q ) 

Alu(K) 

(1/2)A(K) 

(1/2)[K L -B+A'K)] 

(1/2)[2K 0 +A(K)] 

ai l (k) 

-(1/2)A(K) 

(1/2)[K L -B-A'K)] 

(1/2)[2K 0 -A(K)] 

c 

1 

1 

1 

Li (K) J j 

1+(1/2)A(K) 

(1/2)[K L +A 1+ 1+A(K)] 

K L -K 0 +Ai+(1/2)A(K) 

Limiting values for the Kp range for sheared plan forms can be 
determined from the relation: Kp + Kp = 2K 0 

A(K) -- iu(KL-Ku) - 2hi(K L -K 0 ) 


*For K 0 > B the lower limit B is replaced by K 0 . 




APPENDIX D 


SAMPLE CALCULATIONS 


Sample calculations of the dimensionless zero-lift wave drag for 
several of the configurations discussed in this report are presented in 
order to demonstrate the computing procedure in detail. From equa- 
tions (A3) to (A9) the general computing procedure can be stated as 
follows . 

To calculate the drag contribution of each component or the inter- 
ference between each pair of components of an equivalent body of 
revolution 

1. Determine the total length of each component or of each pair of 
components and calculate the moments about the midpoint of this length. 

2. Calculate the length ratios for each component or pair of 
components . 

3* With the information of steps 1 and 2 calculate the quantities 
Lno(P^/ c j 0) ( Gee eq. (A4) ) . 

4. With the information of step 3 it is a simple matter to evaluate 
equations (An) and (A7)j the latter of which represents the dimensionless 
zero-lift wave-drag contribution of one equivalent body. 

9. The results of step 4 for all equivalent bodies are then graph- 
ically averaged according to equation (A9) to obtain the dimensionless 
zero-lift wave drag of each configuration component alone or the inter- 
ference between components, as the case may be. 

If desired, the results of step 4 for the drag of each component and 
the interference between components may be summed before undertaking 
step 5 . Dimensional values are obtained, of course, by evaluating the 
coefficients of equations (All) or (Al6) for airplane configurations or 
systems of bodies of revolution, respectively. 


SINGLE BODY OF REVOLUTION 


From equation (34) the dimensionless area distribution of a Sears- 
Haack body can be written 



The total length is that of the body. Using equation (Dl) in equation (Ala) 
the moments are calculated to oe 


k 

Mok 

k 

Mok i 

0 

3 A 

1 

0 : 

2 

1/y 

pi 

0 ; 

4 

-;/6 4 

5 

0 1 

6 

3/1: :b 

7 

0 ! 


For a single body alone the length ratio is unity. Hence, from equa- 
tion (A17) 

f>20 = ;>A 

anu all other Ln 0 = 0 for n > Therefore, from equations (22) and 
(;].), the dimensionless zero-lift wave drag at Mach number 1 is 

d = u( ;; /i) 2 = 9/a = 3.125 (D2) 


PAIR OF BODIES OF REVOIUTION 


'Hie total zero-lift wave drag will be calculated for the pair of 
Sears -Hauck bodies shown) in sketch (g) for cnly M = 1, since these 
calculations are typical of those for each equivalent body required for 
M > 1. The total zero-lift wave drag for h = 1 is given by (see 
eq. ( 26 )) 


Dp -- Dg + % + Djh ( d 3) 

The first two terms of equation (Df), respectively, represent the drag 
contributions of bodies E and II alone and are given by equation (D2). 

In evaluating the interference term, D 1;II , the dimensionless value 
of the combined length of both bodies is 4. The moments of each body 
about an axis through the origin (see above table for single body of 
revolution) must be transferred to an axis through the midpoint of the 
total length. It can be seen from sketch (g) that such an axis passes 
through the tail of body E and the nose of body H. Hence, in 
dimensionless coordinates, the transfer distances are 


Ae = 1 
A H = -1 


(d4) 



for bodies E and H, respectively. From equation (D4) and from 
the moments calculated above for a single body of revolution as 
of M ok ’ , equation (C19) yields the following moments about the 
of the total length of both bodies . 


use of 
value s 

midpoint 



The ratio of the length of each individual body to the total length of 
bodies is l/m . Hence, from equation (A17) the quantities Ln 0 for both 
bodies are 



Using the above information in equations ( 28 ) and (27) yields 

f>EH = 0.204 (Dp) 

Finally, from use of equations (D ') and (Dp) in equation (D 3 ), the dimen- 
sionless total zero -lift wave drag is 


D*p = 2 . 4p4 


COMPLETE WING 


Calculations of the dimensionless zero-lift wave drag are presented 
for a family of diamond plan-form wings (see sketch (b)) which have 
biconvex airfoil sections of constant thickness ratio. Calculations have 
been made for a Mach number range corresponding to values of 0b/c from 
0 to 3 . 





The important dimensionless quantities describing the plan form, 
the airfoil section, and the spanvlse variation of the thickness along 
lines of constant percent chord are given, respectively, by 

K 0 = A = 0 ; B - (D 6) 

t(cc) = 1 - a 2 (D7) 

$(t)) = 1 - T) (Dei) 

'Plie value of t(a) given by equation (D7) in. actually for a parabolic-arc 
airfoil section. However, for values of the thickness ratio up to about 
o percent, this is a satisfactory approximation . 


Calculations foi* the Wave Drag of Each Wing Panel Alone 


Initially one determines the range of mines of K for which the 
calculations are to he made. Since K 0 ~ 0 and (pb/c) max = the range 
of the values of K for which the moments of one wing panel must be 
calculated are determined from equation (25 ns -• and K.j_ j([1 ^ v - •, 

The moments of each sheared wing panel are first calculated about 
the 1 0 -percent chord line at the plane of s; name try. From equations (Cid) 
Ui.u (DY) 


r 


( 1 -a 2 ) a 1 'da 


% =< 


4 f 1 


K \m+l m+3, 


m e ven 


m odd 


Values of 0 £ m < 7 arc required. 
From equations (Clw) and (D8) 


(D9) 


Ek 


-m+i 


= ^ (I-y)ti 1 


k*-m+r 


di] 


■^k-m+r 


k-m+r+1 k-n+rh 


(DIO) 



Values of 0 <(k-m+r) < 8 are required. From use of the results of 
equation (DIO) in equation (Clp) values of F^ m are determined for 
0 < k < 7 and 0 < in < k. With F^ m and I m determined, moments about 
the £ = 0 axis are determined for the sheared wing panels by use of 
the desired values of “3 £ K < 3 in equation (Cl6). Because the 
airfoil section is symmetrical about the 30 -percent chord, in this case 


Thus , 


(-K) = Mok' (K) for k even 
Mq^ 1 (" K ) = "Mok T (h) for k odd 
only values of 0 < K < 3 are necessary. 


(DU) 


Next, the moment of each sheared 'wing panel alone must be transferred 
to an axis through the midpoint of the total length of each sheared wing 
panel. From the table of appendix C for the sheared wing panels alone, 
the transfer distances and length ratios are 


K range 

A W (K) 

[c/z(K)] W 

K < -1 
-1 < K < 1 
K > 1 

(1/2) (K+l) 
(1/2) (K-l) 

2/(l-K) 

2/(l+K) 


With the above values of Ay(K) and the values of Mop T (K), values of 
Mok(K) are determined from equation (C19) for 0 < k < 7- For this special 
case, the relationships of equation (Dll) also apply to the moments about 
the midpoint of the total length of each sheared wing panel. With the 
information now available the dimensionless zero-lift wave drag of the 
upper and lower sheared wing panels alone can be calculated from equa- 
tions (23), (22), and (21) for 2 < n < 9 . The results are plotted in 
figure 7 a function of K. It will be noted that for this case, 

Dy(K) is symmetrical about K = 0. This is again the result of the 
symmetrical airfoil section. In view of the relationships of equa- 
tion (Dll) it can be seen from equation (A17) that 


Lno(-K) = L no (K) for n even 
L no (-K) = -L no (K) for n odd 


(D12) 


Since [L no (K)] 2 is required (see eq. (22)) and the length ratios are 
symmetrical about K = 0, Dy(K) is necessarily symmetrical about K = 0. 
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Calculations for the Interference Wave Drag Between Wing Panels 

To calculate the interference between wing panels the ordered pairs 
of upper and lower sheared wing panels are first determined from equa- 
tion (jl); thus 

K U = " k L (DIj) 

The moments, Mok'(K), for the ordered pairs of sheared wing panels must 
be transferred to the midpoint of the total length of both sheared wing 
panels. From the table in appendix C for interference between sheared 
wing panels , the transfer distances and length ratios are 


Kl range 

Alu(K) 

ai l (k) 

[c/zooq 

o < Kl < 1 

0 

0 

i ! 

Kl > 1 

0 

i 

0 

i/k l 


From equations (29) and (jO), values of lno(*0 ure calculated for the 
pairs of upper and lower sheared wing panels for 2 < n < 9- Again, as 
a result of the symmetry of the given configuration, use of equation (Dig) 
in equation (D12) indicates that the detailed calculations of values of 
L no (K) for interference are necessary for only the lower sheared wing 
panels. The dimensionless interference zero-lift wave drag between sheared 
wing panels is given by equations (28) anc (27). The results are plotted 
directly as a function of f3b/c cos 0 in figure 8. Also shown in 
figure 8 are the dimensionless zero-lift vave-drag results for each of 
the sheared wing panels alone replotted as a function of pb/c cos 0 
from the data of figure 7* It is emphasi 2 ed that all results shown in 
figures 7 8 were calculated using six significant figures. 

Calculations for the Total Vave Drag of the 
Complete Wing 

The total zero-lift wave drag for each complete sheared configuration 
is calculated using the form of equation (26) for an equivalent body. 

The results are included in figure 8. The variation of the dimensionless 
total zero-lift wave drag of the complete wings ’‘With pb/c is obtained 
by replotting the total wave-drag results of figure 8 against 0, for 
constant values of pb/c, and graphically averaging the results as 
indicated by equation (20). It should be recalled, however, that complete 
plan forms are symmetrical in jt/2. Values of Dp (pb/c, 0) obtained from 
the data of figure 8 for several values o 1 pb/c are shown in figure 9 
replotted against 0. The dimensionless total zero-lift wave -drag results, 
D t ( pb/c ) , are shown in figure 10 . 

The data of figure l(u) were obtained from use of the results of 
figure 10 in equation (33). For this complete plan form b^b = t x /t = 1. 
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APPENDIX E 


EVALUATION OF THE INTEGRAL FUNCTIONS FOR AIRFOIL SECTIONS 
WITH NONANALYTIC THICKNESS DISTRIBUTIONS 


Determination of the moments used in the zero-lift wave-drag 
calculations depends upon the evaluation of the integral functions of 
the airfoil-section thickness distribution given in appendix C by 

In = - J T(a)a m da (C12) 

where t(cx) is the dimensionless airfoil-section thickness distribution 
for -1 < a < 1. Equation (C12) can be evaluated with little or no 
difficulty for airfoil sections for which an analytical expression for 
t(co) is available, such as an NACA four-digit series airfoil section. 

The thickness distribution is not available in analytical form for many 
airfoil sections, such as an NACA six-digit series airfoil section. Tie 
thickness is generally specified, however, at about 20 airfoil chord 
stations. It appears, therefore, that a graphical evaluation of equa- 
tion (C12) is possible. Unfortunately, experience has demonstrated that 
the accuracy of such a technique 
is insufficient to provide reli- 
able wave-drag results. Conse- 
quently, special techniques are 
required to obtain an approximate 
expression for the thickness dis- 
tribution of an arbitrary airfoil 
section, and hence, evaluate equa- 
tion (C12) analytically . 



In appendix C of reference 4 
it was found possible to approxi- 
mate the thickness distribution of 
an arbitrary airfoil section in 
analytical form with a series of 
adjacent parabolic -arc segments of 
continuous slope. This is possible 
when the airfoil-section thickness 
is specified at a sufficient number 
of chord stations, say I, to define 
adequately the thickness distribu- 
tion, and the slope of the thick- 
ness distribution is given for at 
least one chord station (see the 
lower part of sketch (i)). The 
salient feature of this type of 
approximation which is important 
to the evaluation of equation (C12) 


r\a) 





2-0 


is the fact that the second derivative of the thickness distribution is 
constant over each segment of parabolic arc (see the upper part of 
sketch (i)). Hence, one may write for the second derivative of the 
dimensionless thickness distribution 


t" (a) ^ r’^Oi) - constant for 


a i < a < a i+i 


1 < i < (1-1) 


(El) 


The details for finding the analytical expression for t(cc,), t ' (a) , 

and t" (a) can be found in appendix C of reference 4. With this information 

it is possible to evaluate equation (C12) in a simple analytical fashion. 


If one groups the integrand of equation (C12) as t(cx,) and a m da, 
a partial integration yields 


In 


ii/jL 

it |m+l 


[T(l)+(-l) m T(-l)] 


m-tl 



(a)a m+1 da 


■1 


J 


(E2) 


Similarly, if one groups the integrand of the integral expression in 
equation (E2) as r T (a) and a m+1 da, a partial integration yields 


2 1 


nn 


* rn+1 


[T(l) + (-l) m T(-l)] - 


1 


(m+l)(m+2) 


[T'(l) + (-l) m+1 T>(-l)-Pj 


where 


Pm = 


II / \ ITH 2 . rj 

t (aja da 


(E3) 


(e4) 


-1 


As a result of the parabolic-arc approximations of t(<x), t" (a) is 
constant over each segment of arc as indicated by equation (El). Hence, 
with the aid of equation (El), equation (l4) can be written 

I-i 


P m - Y VK) f 1 - 


^ nr. • 


1=1 

which, after integration, becomes 

I-i 

T, ~ 1 


L m 


m+3 


Z tt/ \ f m-3 m+s\ 
T ( a i) (^i-i - a i ) 


(E5) 


1=1 
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(b) Triangular plan form, 










( b) Exposed wing panels of a wing -body combination 



Body E of given and 
sheared configuration 
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Figure 6.- Sketch, in dimensionless coordinates, showing lengths and moment transfer distances 

for sheared exposed plan forms. 



Figure 6.- Continued. 




(C) 0 






Interference between 
sheared wing panels 



Figure 8.- Variation of dimensionless zero-lift wave drag with, (3h/c cos 0 for the sheared 

components of the wings of appendix D. 
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